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I. 


SUMMARY 


A nonlinear  theory  for  solving  partially  cavitating  cascade  problems 
was  formulated  in  two  different  types  of  input  data  set-ups,  one 
specifying  cavitation  number  and  the  other  the  cavity  length.  The 
second  method  was  necessary  to  find  a unique  solution  because  the 
cavity  length  has  two  different  values  for  a specified  cavitation 
number  in  some  blade  geometry  cases.  This  fact,  together  with  the 
investigation  of  detailed  pressure  distributions  on  the  upper  wetted 
part  of  the  blade,  explained  the  unsteady  cavity  phenomena  previ- 
ously reported  in  some  experiments.  Suggestions  were  made  on  how 
to  avoid  the  occurrence  of  the  unsteady  partial  cavity  phenomena 
when  one  designs  fluid  machinery  used  in  a partially  cavitating 
regime . 
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NOMENCLATURE 


scale  factor  for  mapping  between  W and  £ planes. 

^-coordinate  corresponding  to  the  upper  cavity  separation 
point 

^-coordinate  corresponding  to  the  cavity  reattachment  position 
and  also  used  as  a chord  length  of  blade 
lift  and  drag  coefficients  parallel  and  normal  to  the 
incoming  flow  direction 

( = L or  DA  PU^  c ) 

lift  and  drag  coefficients  in  the  direction  of  the 
x and  y axes 

lift  coefficients  in  the  x-direction  normalized  by 
the  geometric  mean  velocity  Um 
pressure  coefficient  |=  | p-pjA  pU^ ) 

drag  force 

spacing  of  cascade  blades 

^-coordinate  corresponding  to  the  upper  trailing  edge 

nonlinear  equations  representing  the  boundary  value 

problem  for  unknown  solution  parameters  x^ 

index  for  complex  variable  | = V-lJ 

lifting  force 

cavity  length 

static  pressure 

flow  velocity  modulus 

total  arc  lengths  of  the  upper  and  lower  wetted  part 
of  the  blade,  respectively 


* 


sg  (x) 


local  arc  length  measured  from  either  the  leading 
edge  or  the  reattachment  point  of  cavity 


sign  function 


= 1 for  x > 0 


=-l  for  x < 0 


thickness  of  plano-convex  blade  normalized  by  chord 
flow  velocity  amplitude 


complex  potential  = 4>+i4'J 
physical  plane  j=  x+iyj 

flow  angle  made  with  the  x-axis  of  the  physical  plane 
local  blade  slope  measured  from  the  x-axis,  clockwise 

negative;  for  the  second  arc  S2,  -tt  is  to  be  added 
geometric  stagger  angle 

stagger  angle  of  the  potential  plane  cascade  set-up 

(=  a1+Y] 

transform  potential  plan  (=  £+iri| 


^-coordinate  corresponding  to  the  upstream  infinity 

of  the  physical  plane 

vertical  coordinate  in  the  £-plane 

flow  angle 

horizontal  coordinate  in  the  £-plane 


cavitation  number 


(pl-pc)APuI 


logarithmic  velocity  normalized  by  U2 
potential  function 
stream  function 


= £n  (q/U2 )■ 


hodograph  variable  j = 9+ix 
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INTRODUCTION 


It  is  a common  practice  that  hydrodynamists  avoid  high  speed 
fluid  machinery  or  components  to  be  used  under  partially  cavi- 
tating  conditions  because  many  adverse  hydrodynamic  effects  exist 
for  such  conditions.  One  of  the  most  important  phenomena  is  that 
of  cavitation  damage  caused  by  ca /itv  bubble  collapse  taking 
place  near  the  material.  This  material  damage  is  so  severe  that 
components  such  as  marine  propellers  after  several  hundred  hours 
of  exposure  to  cavity  collapse  will  become  totally  out  of  use 
(see  e.g.  [l]  for  such  examples).  Furthermore,  it  was  found 
in  experiments  that  partially  cavitating  flows  experience  inher- 
ent unsteady  flow  phenomena;  the  Length  of  cavity  changes  from 
a very  short  one  to  one  chord  ler  jth  in  an  oscillatory  way  (see 
e.g.  [2]  and  [3]).  It  is  for  this  reason  that  effort  by  the 
designer  is  directed  toward  either  totally  avoiding  cavitating 
conditions,  if  this  is  possible,  or  purposely  generating  super- 
cavitating  conditions  in  which  cavity  length  is  larger  than  the 
body  and  thus  bubble  collapse  occurs  behind  the  material. 

However,  in  recent  years  some  marine  and  fluid  machinery  have 
quite  often  been  used  in  the  partially  cavitating  flow  regime. 

One  such  example  is  in  the  marine  propeller  area;  convential  and 
supercavitating  propellers  are  sometimes  used  under  off-design 
conditions,  the  former  at  a lower  advance  speed  and  the  latter 

* Numbers  in  brackets  designate  references  at  end  of  paper. 


at  a higher  advance  speed  than  that  of  the  design  point.  In 
both  cases  some  portions  of  propeller  blaues  have  partial  cavi- 
ties as  were  shown  in  experimental  works  by  Sontvedt  [4],  Bohn 
[5]  and  Peck  [6]. 

There  exists  more  positive  designers 1 attitude  of  using  the  par- 
tially cavitating  condition  in  some  area  of  fluid  m’chinery  such 
as  pump  inducers.  Cavitating  inducers  used,  for  example,  in 
liquid  propellant  rockets  always  operate  at  a partially  cavitat- 
ing condition.  Because  of  simplicity,  such  inducer  blades  are 
made  of  helical  profile  having  long  straight-line  face.  The 
cavity  length  is  designed  to  be  always  shorter  than  the  chord 
length  (see  the  work  of  Acosta  [7]). 

Interest  has  therefore  grown  recently  in  how  to  analytically 
determine  the  force  coefficients  of  partially  cavirating  single 
foil  or  cascade  blades.  Linearized  theories  were  applied  by 
many  researchers  for  both  single  and  cascade  blades.  The  former 
include  works  by  Acosta  [8],  Geurst  [9],  Geurst  and  Verbrugh 
[10],  Wade  [ll]  and  Golden  [12].  Acosta  [8]  and  Geurst  [9] 
treated  the  flat  plate  single  foil  problem  and  Geurst  and  Verbrugh 
[10]  solved  the  problem  for  camber  line  profiles.  Wade  [ll]  and 
Golden  [12]  considered  both  camber  and  thickness  effects  in  the 
same  linearized  theory  framework  but  with  different  solution 
methods;  the  former  used  an  analytical  solution  method  for  vor- 
tex and  source  distributions  whereas  the  latter  employed  a 
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numerical  method  to  determine  the  strength  of  the  singularities. 
For  cascade  configurations,  Wade  [13 ] solved  flat  plate  cascade 
problems  in  exactly  the  same  manner  as  for  single  foils.  Algebra 
for  the  cascade  flow  becomes  more  complicated  due  to  extra  cas- 
cade mapping  necessary  for  this  case. 


As  has  been  seen  above,  all  the  existing  analytical  methods 
regarding  the  partially  cavitating  flow  employed  the  linearized 


theories.  The  limitations  of  such  theories  are  well  described 
in  each  paper  [8-13];  any  perturbed  flow  quantities  due  to  the 
existence  of  cavity  and  obstacle  in  the  uniform  flow  should  be 
small.  Ranges  of  flow  incidence  angle,  cavitation  number  and 
camber  of  blades  for  which  the  accuracy  of  the  linearized 
theory  is  maintained  to  a desirable  degree  are  not  known  until 
comparisons  with  experiments  are  made.  Questions  particularly 
arise  for  the  linearization  procedure  for  the  nonlinear  boundary 
condition  on  the  cavity  wall.  Unlike  the  supercavitating  flow 
regime,  the  cavitation  number  o is  relatively  large  compared  to 
the  uniform  flow  velocity.  Nevertheless,  all  the  linearized 
theories  mentioned  above  assumed  a to  be  small.  Furthermore, 
partially  cavitating  flow  around  the  cascade  having  relatively  high 
solidity  will  generate  large  disturbed  velocity  components  for 
which  the  application  of  the  linearized  theory  is  no  longer 
valid. 


Under  the  circumstances  just  described,  an  analytical  tool  for 
determining  the  hydrodynamic  forces  of  partially  cavitating 
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cascade  blades  has  been  developed.  The  method  used  here  is 
the  nonlinear  potential  theory,  covering  a complete  range  of  geo- 
metric configurations  and  flow  parameters  without  generating  any 
inaccuracies  in  satisfying  the  boundary  conditions.  The  method 
is  quite  similar  to  that  used  for  supercavitating  cascade 
flow  analysis  by  Furuya  [14],  but  with  substantial  increase  in 
complexity  both  in  mathematical  formulation  and  numerical  com- 
putations . 

In  the  following  sections,  mathematical  formulations  and  numer- 
ical methods  will  be  described.  Some  representative  numerical 
results  will  be  compared  with  those  of  a linearized  theory  [13] 
and  an  experiment  [3]. 


2. 


NONLINEAR  BOUNDARY  VALUE  PROBLEM 


A mathematical  model  representing  the  physical  flow  of  partial 
cavity  around  the  cascade  has  to  be  established  prior  to  the 
formulation  of  a nonlinear  boundary  value  problem . The 
reattachment  of  cavity  on  the  upper  portion  of  the  blade  makes 
such  a flow  model  construction  difficult.  Highly  turbulent 
flows  near  the  end  point  of  cavity  cannot  be  accurately  repre- 
sented within  the  framework  of  potential  flow  assumption  cur- 
rently taken.  Among  many  possible  flow  models,  we  have  chosen 
a most  realistic  one  as  shown  in  Figure  1.  The  cavity  streamline 
detached  from  the  leading  edge  of  the  blade  terminates  at  point 
C,  the  location  of  which  is  either  specified  or  determined  as 
part  of  the  solution  as  a function  of  cavitation  number.  This 
depends  upon  the  type  of  input  for  boundary  conditions,  i.e. , whether 
H or  a is  used  as  input  data.  Then  the  same  streamline  follows 
the  body  profile  $>2  (x)  through  the  trailing  edge.  The  gap  existing 
between  the  streamline  and  the  upper  body  boundary  is  considered 
to  be  a turbulent  wake  region,  which  is  implicitly  eliminated  from 
the  present  potential  theory  analysis.  From  the  trailing  edge, 
the  wake  with  a constant  pressure  condition  extends  to  the  down- 
stream infinity  and  is  assumed  to  close  there. 

The  end  point  of  cavity  is  interesting  from  both  physical  and 
mathematical  viewpoints.  In  the  present  flow  model,  the  point 
C is  termed  as  a single  spiral  vortex  model  named  after  Tulin 
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[16].  From  the  cavity  side  the  streamline  curls  in  to  a point  C, 
then  immediately  follows  the  corresponding  body  slope.  It  seems 
that  this  behavior  of  the  streamline  more  or  less  represents  the 
cavity  collapse  in  the  actual  flow  situation. 

It  must  be  noted  that  this  type  of  cavity  closure  model  may  not 
be  suitable  for  very  thick  foil  cases,  since  the  gap  between  the 
streamline  and  body  boundary  mentioned  above  is  out  of  control. 

It  may  happen  that  the  cavity  ends  inside  the  blade  structure 
if  the  blade  thickness  is  large.  For  such  cases,  one  has  to  use 
an  open  wake  cavity  model  in  which  the  cavity  closure  condition 
at  down  infinity  is  relaxed  but  instead  replaced  by  specifying 
the  cavity  thickness  at  the  end  point  of  cavity.  In  this  way 
the  streamline  crossing  problem  mentioned  above  will  be  totally 
eliminated.  The  mathematical  formulation  for  the  open  wake  model 
is  no  more  complicated  than  the  closed  wake.  In  this  paper,  how- 
ever, the  latter  will  be  used  throughout  investigation. 

The  mathematical  formulation  is  quite  similar  to  that  of  the 
supercavitat ing  cascade  flow  problem  in  [14],  thus  repetition 
of  detailed  explanations  will  be  avoided. 

The  partially  cavitating  flow  around  the  cascade  shown  in 
Figure  1 is  first  mapped  onto  the  potential  plane  W = <f+i'i'  and 
then  to  the  upper  half  of  the  transform  plane  c by  the  cascade 
mapping  function,  see  Figures  2(a)  and  (b) , respectively.  The 
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cascade  mapping  function  and  its  derivative  are  given  by  the 
following  equations: 


w - & 


i5 , 


(1-5/?^  + e inll-C/^  ) 


(1) 


where 


C.  = A e 


i (tt/2-6  ) 


<5  = o^+y. 


(2) 

(3) 


and  A is  a scaling  factor  between  the  two  planes,  yet  to 
be  known.  Thus,  the  derivative  of  W on  the  5 axis  is  given 

5 cos  5 


dW 

d? 


d 

ir 


(5-AsinS)^  + (A  cos  5)^. 


(4) 


The  stagger  angle  6 of  the  W-plane  cascade  set-up  in  Equation  (3) 
was  determined  by  the  coordinate  point  (<*>g,  >Fq)  where 
*0  = du^*  sin  (a1+y) 

'•'g  = dU1’cos  (a1+y)  . 

Introducing  a hodograph  variable  id 

m = q e-i9  = U,e-iu 
dz  2 


(5) 


(6) 


or 

(D  » 9 + i T 
r - ln(q/U2) 

where  U2  denotes  the  velocity  at  downstream  infinity,  we  can 
express  the  boundary  conditions  on  the  real  axis  £ in  the  Q -plane 
as  follows: 


(7) 

(8) 
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(i) 


t = 0 


< S < - 1 


f < S 


< 00 


t = in  (VI  + a/U2)  / b < £ < c 


(iii)  9 = 31(=tan  dy/dx) 


- 1 < £ < 0 


3-,  +TT 


32+^ 


0 < £ < b 


c < £ < f. 


This  is  a typical  mixed-type  boundary  value  problem,  the  solu- 
tion of  which  can  readily  be  found  by  the  same  approach  as  used 
in  [14],  but  in  a somewhat  more  complicated  form  due  to  the 
existence  of  two  branch  cuts,  - 1 < £ < b and  c < £ < f. 


Analytically  continuing  the  upper  half  of  £ -plane  into  the  lower 

* * 

half  plane  with  <o(£  ) = co  (£)  , where  * denotes  the  complex  con- 
jugate, these  boundary  conditions  (i)  - (iii)  will  be  expressed 
as  follows: 


(jJ  - (Jj 


=0;-®  <£<-1,  f < £ < 


(ii)  co  - <o  = 


= 2i!n^Vl+a/U2^;  b 


< £ < C 


(iii)  u + 


23-,  ;-l<£<0 

2 (81+ir)  ,•  0 < £ < b 

2 ( 62+1T ) ; c < £ < f . 


The  homogeneous  part  of  the  solution  H(£)  for  which  the  right 
hand  sides  of  (i)  - (iii)  are  all  zero  is  obtained  by  consider- 
ing the  singularity  at  the  end  point  of  cavity  and.  no  singular- 
ities anywhere  else; 
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HU) 


(9) 


U-t-1)  U-b)  U-f) 
c-c 


The  general  solution  G(c)  satisfying  the  mixed  boundary  conditions 
is  defined  by  GU)  = wU)/HU)  and  obtained  by  applying  the 
Cauchy  integral  formula; 


G(C) 


Therefore , 
u>U)  = H(c)G(c) 


2.1  BOUNDARY  VALUE  PROBLEM  TYPE  I (SPECIFY  CAVITATION  NUMBER  o) 


At  the  upstream  infinity, 
<*>  ( C i ) = o1  + i tn  (1/U2)  : 


2 equations 


(ii)  At  the  downstream  infinity, 
id  («)  = a2  : 


1 eauation 


(iii)  Length  of  First  Arc  = S^: 


(iv)  End  point  of  the  streamline 
(£=f)  matches  the  physical 
trailing  edge: 

(v)  Continuity  Equation: 


1 ecuation 


1 equation 
1 equation 


where  and  S2  are  depicted  in  Figure  2 (c) . It  is  noted  that 
the  arc  length  condition  (iv)  replaces  the  potential  difference 
condition  in  the  supercavitating  cascade  prol 'em.  (see  Refer- 
ence [14]) 


Application  of  these  conditions  to  Equation  (10)  yields  the 
following  six  equations; 

i)  At  Upstream  Infinity: 

i 

j'1}  3 fU1+l)U1-b)U~|  k (b 

* ^ cl“c  * J I (s’+l)  (b-e  ■)  U’-fT 


(C+l)  (b-g1)  (S'-f) 
5 1 -c 


d5' 

r-c- 


1 dC 

U’+l)  (b-C  1 ) (g'-f)  5’~ 

C'-c 


+ i ,n  / Vl>c\ 

-?1  * \ U2  / 


J J (S'+l)  (g'-b)  (g  *-f ) 

b T 


dg'  . 1 f 32+7r 

'"Cl  ' l J (g’+l)  U~ 


b) (f-51)  " ’1 


C ' -c 


_ ) a^+i  in(l/U2)|  = 


u>(0 


-1  < 5 < 0 
0 < £ < b 


The  arc  length  condition  is  therefore  satisfied  by  the  follow- 


ing equation; 

f4=  s!  ' = 0 (14) 

where  dW/d£  is  given  by  Equation  (4). 


iv)  End  Point  of  the  Streamline  at  g = f Matching  the  Trailing  Edge 
of  the  Blade: 


Similar  to  the  above  derivation,  the  equation  f is  obtained  as 
follows : 


(f)  = 0 


(15) 


f5  H S2  s 2 


where 

s. 


(5)  - f €liU)  |W  as. 

J u2  dg'  ’ 


and 
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g2UJ  ^ 


-*/ 
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iJ 

-/ 


J (g’+l)  (b-C')  (g’-f) 
" 5 '-c 
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dg' 


dg  ' 

g '-g 


■ Hr ) 


f 

1 

dg1  if 

J 

J (5’+l) (g'-b) (g’-f) 

g'-g  kJ 

b 

Y g >-c 

c 

S 2 ■*“Tr 


(g'+l) (g ' -b) (f-g1) 
g ’-c 


f 

V 


dg' 

g'-g 


•(16) 


In  order  to  determine  the  second  arc  length  S2 , the  end  point  of 

cavity  should  be  known.  The  cavity  shape  is  expressed  in  terms  of 
g -coordinate; 


x-x 


B 


y-yB 


b 


gc(5f) 


9c(C') 


dW  dg’ 
dg' 


(17) 


(18) 
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where 


gcu) 


(g+1)  (C~b)  (f- 

c-g 


B,  (g  ' ) 


+1)  (b-C ' ) (g'-f) 


dg  ' 
g '-g 


g'-c 


-f  ===== 

/ / (g'+i)  (b-e ' 

0\  r ' -c 


dg'  ,1  ln/\Jl+a  \ „ 

) (g'-f)  5'-C  if  \ u2  / 


< 

/ 


vn: 


1 dg ' 

+1)  (g'-b)  (g’-f)  g'-g 


-*/ 


62+u 


g'-c 


dgj_ 

+1)  (g'-b)  (f-g ' ) 


g ' ~c 


x_,  yn  are  the  physical  coordinates  of  the  blade  leading  edge. 

The  end  point  of  the  cavity  (xc,  yc)  is  obtained  by  setting  the 

upper  bound  of  the  above  integrals  (17)  and  (18)  to  be  c.  Once 

(x  , y ) are  determined,  the  second  arc  length  S_  is  easily 
c c ^ 

calculated  by  following  the  body  profile. 


v)  Continuity  Equation: 


(19) 


The  mass  flow  condition  for  the  upstream  and  downstream  is  given 
as ; 


f6=  cos(a^+Y)*d  - CJ2  cos(a2+y)*  d =*=  0 


(20) 


Equations  (11),  (12),  (14),  (15)  and  (20)  now  provide  six 

independent  relationships  for  six  unknown  parameters. 
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2.2  BOUNDARY  VALUE  PROBLEM  TYPE  II  (SPECIFY  CAVITY  LENGTH  Z \ 

— - . — sJL 

The  cavity  length  l is  used  as  input  data  in  this  case  and 
the  corresponding  cavitation  number  a is  determined  as  a result 
of  the  solution.  Minor  changes  to  the  TYPE  I problem  are  neces- 
sary. In  addition  to  in  Equation  (15) , one  more  equation, 
fj,  is  added  in  order  to  satisfy  the  fixed  cavity  length  con- 
dition. 

f7  5 xo  - ‘o  * 0 <21) 

where  xc  is  calculated  by  setting  £ to  be  c in  Equation  (17) . 

We  have  now  seven  independent  equations  for  seven  unknown  solu- 
tion parameters  including  a. 

There  is  absolutely  no  difference  in  final  solutions  between 
the  above  two  boundary  value  problems  as  long  as  the  physical 
and  geometric  flow  conditions  are  the  same.  In  mai.y  actual  cases, 
the  cavitation  number,  which  can  be  simply  determined  by  the 
upstream  static  pressure  and  vapor  pressure,  is  a specified  ijj- 
put  rather  than  the  cavity  length.  Only  the  first  type  of 
B.V.P.  is  therefore  usually  used  to  provide  solutions 
for  the  problem  like  supercavitating  flows.  Special  reasons 
for  introducing  the  second  type  of  B.V.P.  in  this  study  exist: 
i)  Increase  the  stability  of  the  numerical  iterative 
method  by  fixing  the  wetted  portion  of  the  upper 
blade,  and 

ii)  Avoid  the  numerical  instability  caused  by  the  existence 
of  multiple  cavity-length  solutions  for  one  cavita- 
tion number  as  reported  in  [9],  [13]  and  [15]. 
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As  will  be  discussed  in  the  following  section,  the  numerical 
instability  actually  occurred  for  certain  cases  when  the  first 

type  of  B.V.P.  was  used,  whereas  cavitation  number  was  uniquely 
determined  for  any  specified  cavity  length  when  the  second 
method  was  employed. 
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3 . NUMERICAL  METHOD 

The  same  Newton's  functional  iterative  method  as  that  in  [14] 
was  applied  to  the  present  system  of  nonlinear  functional 
equations  for  determining  the  solution  parameters.  Representing 
these  equations  and  unknown  parameters  as 


f (x)  = 0,  ( 

the  iterative  method  is  introduced  in  the  following  relation; 


(x  (n))  • (x(n+1)-  x(n))=  -f(x(n)) 


where  J is  a partial  derivative  matrix  whose  component  is 
expressed 


By  assuming  a starting  set  of  x^  values,  x(n+1)  will  be 
calculated  based  on  Equation  (23) . As  has  been  seen  in  the  preceding 
section,  f also  requires  the  information  regarding  3^  and  as 
a function  of  £.  Although  the  body  inclinations  3^  and  are 
specified  in  the  physical  plane,  we  do  not  know  them  as  a function 
of  mapped  potential  coordinate  £ until  the  complete  problem  is 
solved.  This  inherently  nonlinear  implicit  problem  is  conveniently 
resolved  as  a part  of  the  above  iterative  method  by  gradually 
updating  0^  and  6 2 • The  first  assumption  for  0^  and  &2  is  naturally 
that  of  constant  values,  i.e.,  straight  line  body  assumption. 

These  values  are  then  updated  whenever  the  new  set  of  x values  are 
obtained.  The  iteration  procedure  will  be  continued  until  x as  well 
as  the  functional  relations  for  0^  and  02  against  £ converge. 
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4.  PRESSURE  DISTRIBUTIONS,  FORCE  COEFFICIENTS 

Once  a convergent  solution  is  obtained,  the  calculations  for  Cp, 

CT  and  C_.  are  straightforward. 

L>  D 


Since 


S “ 1 


- (</<s) 


where 


2 _ 2 -i  (w-o)  ) 

q = U2e  » 

0^  = 1 - U 2 exp  Im  L(C) 


ri<e: 


for  - 1<£<  b 


£n  Vi+o 

U-, 


b <E,<  c 


g2U) 


c<£<  f 


The  functions  g (?)  and  g2  (C)  were  already  introduced  in  Equations 
(14)  and  (15)  . 

The  lift  and  drag  coefficients  in  the  direction  of  the  x-  and 


y-  axes  are  given 


r~  - . f c — — 

L ~ J Lp  ds  ' dC 


CD  = 
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where 


dx 

ds 


dy 

ds 


ds 
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r dx 


ds. 


dx 
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for  the  second  arc  S. 
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first  arc  S. 
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for  the 
second  arc  S. 


(c<C<f] 


The  lift  and  drag  coefficients  CL  and  CD  in  parallel  and  normal  to 
the  incoming  flow  angle  are  thus  calculated 
CL  = cos  a1  - sin 

CD  = sin  a1  + cos 


(27) 

(28) 
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5-  RESULTS  AND  CONCLUSIONS 

The  present  theory  was  first  applied  to  partially  cavitating 
flat-plate  cascades  and  the  results  were  compared  with  those  of 
a linearized  theory  of  Wade  [13]  where  the  same  parameters  in 
[13]  are  used  in  Figure  3.  Not  only  the  lifting  force  itself 
but  also  the  trend  of  the  curves  are  significantly  different 
between  the  two  theories.  The  discrepancies  particularly  become 
larger  as  the  cavitation  number  increases.  As  has  been  mentioned 
before,  it  is  considered  that  these  discrepancies  are  attributable 
to  the  loss  of  the  linearized  theory's  accuracy  for  larger  cavi- 
tation number  flows.  Furthermore,  the  incidence  angle  taken  in 
the  nonlinear  theory,  i.e.,  8 degrees,  may  be  already  beyond 
the  applicable  range  of  the  linearized  theory. 

The  theory  was  also  compared  with  some  experimental  data  of  par- 
tially cavitating  cascade  of  plano-convex  blades  [3].  The  blade 
thickness  used  for  the  experiment  was  8 percent  chord.  In  order 
to  gradually  update  the  present  numerical  calculations,  the 
thickness  of  plano-convex  blade  was  changed,  starting  from  0 
percent  to  2,  5 and  8 percents.  Each  computation  utilized  con- 
verged solution  of  previous  results  as  a starting  point  of  the 
iterative  procedure.  The  lift  and  drag  coefficients  are  shown 
in  Figures  4 and  5.  Convergent  solutions  were  obtained  for  0 
and  2 percent  thickness  cases  with  the  TYPE  I B.V.P.  in  which 
the  cavitation  number  was  specified.  The  results  not  only 
check  with  the  lift  coefficients  of  fully  wetted  flat  plate 
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cascade  calculated  based  on  Thwaites  [17],  but  also  show  the 
effects  of  blade  camber  as  an  increase  in  The  calculated 

results  of  2 percent  thickness  approach  closer  to  those  of  the 
experiment  having  blades  of  8 percent  thickness. 

We  encountered  severe  numerical  instability  with  the  same  computer 
program  when  the  blade  thickness  was  increased  to  5 and  8 percents. 

No  matter  how  carefully  initial  values  were  chosen,  the  length 
of  cavity  at  each  iteration  wildly  changed  over  a total  range  of 
chord  length  and  never  converged.  It  was  this  time  that  the 
TYPE  II  B.V.P.  approach  in  which  the  cavity  length  is  specified 
and  the  corresponding  cavitation  number  is  determined  was  used. 

The  problem  of  numerical  instability  was  resolved  with  this  second 
method.  As  was  suspected  from  the  occurrence  of  numerical  in- 
stability for  thickness  of  5 percents  with  the  first  method,  it 
was  found  here  that  the  cavity  has  two  different  lengths  for 
one  cavitation  number  as  shown  in  Figure  6.  The  lift  and  drag 
coefficients  obtained  are  also  plotted  in  Figures  4 and  5 and 
showed  peculiar  behavior  in  the  - a and  CD  - o curves.  As 
the  cavity  length  decreases,  the  cavitation  number  first  increases 
but  then  starts  decreasing  with  a slight  increase  of  C^.  It 
seems  that  these  multiple  cavity  length  solutions  and  peculiar 
force-vs.-o  behaviors  will  well  explain  the  unsteady  oscilla- 
tory phenomena  of  the  partially  cavitating  single  foil  or  cas- 
cade flows  as  were  observed  in  experiements  [2]  and  [3].  A 
recent  study  of  Uhlman  and  Jiang  [15]  attributed  the  existence 
of  two  cavity  lengths  for  a specified  a to  the  inaccuracy  of 
the  linearized  theory  and  the  inherent  flow  instability.  The 
first  cause  is  now  obviously  denied. 

J 
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The  plano-convex  foils  such  as  of  0 and  2 percent  thickness 


have  a single  value  of  cavity  length  for  a specified  o while  as 
the  thickness  increases  to  e.g.,  5 percents  the  cavity  length 
starts  having  double  values.  This  fact  seems  to  provide  an 
essential  key  for  understanding  the  above  unsteady  cavity  phenomena. 
The  detailed  pressure  distributions  on  the  upper  blade  surface 
are  plotted  in  Figures  7 (a)  to  (c) . For  the  thickness  of  0 and 
2 percent  cases,  it  is  seen  from  these  figures  that  the  pressure 
distributions  on  the  upper  solid  boundary  always  indicate  positive 
values  relative  to  the  cavity  pressures  Cpc.  On  the  other  hand, 
for  the  blade  thickness  of  5 percent  (Figure  7 (c))  , these  pressures 
rapidly  become  negative  relative  to  CpC,  particularly  as  the  cavity 
length  becomes  shorter.  This  violates  a theorem  applied  to  the 
present  boundary  value  problem  that  the  minimum  pressure  exists 
inside  the  cavity  (see  e.g.,  a textbook  by  Birkhoff  and  Zarantonello 
[18])  . Furthermore,  even  taking  viscous  effects  into  consideration, 
high  negative  pressures  cannot  be  sustained  in  actual  flow  field, 
indicating  that  the  cavity  is  ready  to  extend  to  these  lower 
pressure  regions.  As  soon  as  the  cavity  becomes  longer,  no  more 
low  pressure  field  exists.  As  is  seen  in  Figure  7 (c) , e.g., 
o = .37  case,  the  vaporization  process  is  not  fast  enough  to  pro- 
vide vapor  to  hold  the  cavity.  The  cavity  bubble  is  thus  shed 
downstream  and  a new  short  cavity  starts  growing. 

The  generation  of  higher  negative  pressures  for  thicker  blades 
is  easily  understood;  the  higher  the  negative  curvature  of  the 
obstacle  in  the  flow  field,  the  lower  the  negative  pressure. 

As  a matter  of  fact,  the  present  explanation  for  the  unsteady 
partial  cavity  flow  motion  based  on  the  steady  state  flow  analysis 


22 


is  consistent  with  experimental  observations.  Experiments  of 
Meijer  [2]  and  Wade  and  Acosta  [3]  are  two  well-known  cases  in 
which  the  unsteady  phenomena  were  observed.  Although  one  tested  a 
single  foil  and  the  others  cascade  foils,  a most  significant  coin- 
cidence between  the  two  is  that  both  used  blades  whose  upper  profile 
is  of  unfavorable  shape  for  the  pressure  field,  i.e.,  concave  shape, 
in  these  cases,  circular  arc.  It  is  now  clearly  explainable  why 
the  unsteady  cavity  oscillations  occurred  in  both  cases. 
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Based  on  the  above  discussions,  it  is  realized  that  the  partially 
cavitating  flow  is  either  stable  or  unstable,  depending  upon  the 
pressure  distributions  on  the  we tted  part  of  the  upper  blade  sur- 
face. Figure  6 indicates  the  stable  partial  cavity  condition 
for  thin  flat-plate  cascades.  In  the  past,  little  hydrodynamic 
considerations  have  been  given  for  design  of  blade  backside  pro- 
file. For  example,  the  backside  profiles  of  supercavitating 
propellers  were  determined  from  a structural  integrity  point  of 
view  within  the  limit  of  cavity  envelope.  In  order  to  avoid  the 
occurrence  of  cavity  oscillations  if  such  propellers  operate  at 
an  off -design  point  in  the  partially  cavitating  regime,  it  is 
recommended  that  the  blade  backside  profile  be  properly  designed 
so  that  favorable  pressure  distributions  Cp  prevail.  At  least, 

Cp  on  the  backside  of  blades  under  several  partially  cavitating 
conditions  must  be  investigated  during  the  design  procedure;  it 

should  be  ensured  that  C is  much  larger  than  -a  over  the  full 

P 

range  of  wetted  upper  boundary. 


j 
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Unlike  linearized  partial  cavity  or  super  cavitating  theories, 
the  present  nonlinear  partial  cavity  theory  has  no  singular  be- 
haviors as  Zc  approaches  the  chord  length  nor  does  the  nonlinear 
super-cavity  theory  [14].  In  Figures  4 and  5,  the  regimes  as  ic 
becomes  c are  depicted  by  dashed  lines.  This  does  not  mean  that 

and  CD  are  unknown  or  become  infinite,  but  they  simply  have  slight 
discrepancies  between  the  two  nonlinear  theories  because  the  flow 
models  used  for  the  partial  and  super-cavity  regions  cannot  be 
smoothly  blended. 

Figure  8 shows  the  boundary  profiles  of  typical  partial  cavities 
for  = 9°,  y = 36°  and  blade  thickness  of  2 percents  at  various 
cavitation  numbers  or  cavity  lengths. 

Finally,  it  must  be  noted  that  the  numerical  iterative  method 
used  here  provided  stable  convergent  solutions  as  long  as  one 
chose  a right  type  of  input  data  set-up  as  has  been  mentioned 
above.  Computer  execution  time,  however,  was  found  to  be  two 
to  six  times  that  of  supercavitating  cascade  computations  in 
[14]. 
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Figure  7(c)  The  same  as  Figure  7(a)  excep 


Profiles  of  calculated  cavity  boundaries  for  cascade  of  plano-convex 
blades  having  2%  thickness  at  n1=9°  with  y=36°  and  solidity  = .625 


